Kinetic-growth self-avoiding walks on small-world networks 
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Kinetically-grown self-avoiding walks have been studied on Watts-Strogatz small-world networks, 
rewired from a two-dimensional square lattice. The maximum length L of this kind of walks is 
limited in regular lattices by an attrition effect, which gives finite values for its mean value (L). 
For random networks, this mean attrition length (L) scales as a power of the network size, and 
diverges in the thermodynamic limit (system size N — » oo). For small- world networks, we find a 
behavior that interpolates between those corresponding to regular lattices and randon networks, for 
rewiring probability p ranging from to 1. For p < 1, the mean self-intersection and attrition length 
of kinetically-grown walks are finite. For p — 1, (L) grows with system size as N 1 ' 2 , diverging 
in the thermodynamic limit. In this limit and close to p = 1, the mean attrition length diverges 
as (1 — p)~ 4 - Results of approximate probabilistic calculations agree well with those derived from 
numerical simulations. 

PACS numbers: 89. 75. He Networks and genealogical trees - 05.40.Fb Random walks and Levy flights - 
89.75.Fb Structures and organization in complex systems - 87.23.Ge Dynamics of social systems 



I. INTRODUCTION 

In the last few years, researchers have been accumulat- 
ing evidence that several kinds of real-life complex sys- 
tems can be described in terms of networks or graphs, 
where nodes represent typical system units and edges 
represent interactions between connected pairs of units. 
This topological description has been applied to study 
both natural and man-made systems, and is currenlty 
employed to describe various types of processes taking 
place on real systems (social, biological, economic, tech- 
nological) d El IE 0, M ■ Two highlights of these devel- 
opments were the Watts-Strogatz small-world model [6[ 
and the so-called scale- free networks p|, which incorpo- 
rate several basic ingredients of real systems. 

In particular, Watts-Strogatz small-world networks are 
well suited to study properties of systems with underly- 
ing topological structure ranging from regular lattices to 
random graphs [a|9(, by changing a single parameter [ig |. 
These networks are based on a regular lattice, in which a 
fraction p of the links between nearest-neighbor sites are 
randomly replaced by new random connections, thus cre- 
ating long-range "shortcuts" and causing an important 
decrease in the average distance between nodes @, E| . In 
small-world networks one finds at the same time a local 
neighborhood (as for regular ^lattices) and some global 
properties of random graphs [8[ . In particular, the "dis- 
tance" between any two elements is small as compared 
with the system size, and the propagation of information 
(or signal, disease, damage, .J takes place much faster 
than in regular lattices 0, El, El EE El EE EE El- 
This short global length scale affects strongly the behav- 
ior of statistical physical problems on small-world net- 
works, such as site and bond percolation [18l Qj| > as well 
as the Ising model [H El El- 

Social networks form the substrate where several dy- 
namical processes take place. These networks have the 



property of being searchable, i.e. people (nodes) can di- 
rect messages through their acquaintances to reach dis- 
tant specific targets in only a few steps [23|, Ej, Efl • Al- 
though some dynamical processes on complex networks 
have been studied by using stochastic dynamics and ran- 
dom walks [2y, Ell EM] > it is known that real navigation 
and exploration processes are usually neither purely ran- 
dom nor totally deterministic 0, EiJ, ISO]. In this con- 
text, the generic properties of deterministic navigation 
[31| and directed self-avoiding walks J29[ in complex net- 
works have been analyzed in the last years. Self-avoiding 
walks (SAWs) may be more effective than unrestricted 
random walks in exploring a network, as they cannot re- 
turn to sites already visited. This property has been used 
to define local search strategies in scale- free networks [3fJ] . 
Nevertheless, the self-avoiding property causes attrition 
of the paths, since a large fraction of paths generated 
in a stochastic way have to be abandoned because they 
overlap themselves. This can be an important limitation 
for exploring networks with this kind of walks. 

SAWs have been employed for many years for modeling 
structural and dynamical properties of macromolecules 
[33 . l33j , as well as for studying critical phenomena in lat- 
tice models [3J, |35[ . They are also useful to characterize 
complex crystal structures [36| and networks in general 
[371 ] . In particular, the asymptotic properties of SAWs on 
small- world networks were studied in connection with the 
so-called connective constant or long-distance effective 
connectivity [38|. Recently, kinetic-growth self-avoiding 
walks on uncorrelated complex networks were considered, 
with particular emphasis upon the influence of attrition 
on the maximum length of the paths [39|, |4f|. It was 
found that the average length scales as a power of the 
system size N, with an exponent that depends on the 
characteristics of the considered networks. For regular 
lattices, however, the average maximum length of kinetic 
growth SAWs is finite H [42J , contrary to uncorrelated 



networks, for which it diverges in the thermodynamic 



limit (N — > oo). Thus, for small- world networks there 
will appear a crossover when changing the rewiring prob- 
ability p, from a regime with confined paths (finite av- 
erage length typical of regular lattices, p = 0) to one 
with diverging length (characteristic of random networks, 

p-»l). ' ' ' 

Here we study long-range properties of kinetically- 
grown walks on small-world networks, and consider the 
"attrition problem" . The number of surviving walks to 
a given length n is obtained by approximate analytical 
procedures, and the results are compared with those ob- 
tained from numerical simulations. We find that the aver- 
age maximum length increases as the rewiring probability 
p is raised, and eventually diverges in the thermodynamic 
limit for p — > 1, as expected for random networks. We 
note that the term length is employed throughout this 
paper to indicate the (dimensionless) number of steps of 
a walk, as usually done in the literature on networks [3| , 
and does not correspond to a distance on the underlying 
regular lattice. 

The paper is organized as follows. In Sec. II we 
give some definitions and concepts related to SAWs, 
along with details on our numerical method. In Sec. Ill 
we analyze the length at which non-reversal random 
walks intersect themselves on small-world networks (self- 
intersection length), and in Sec. IV we calculate the av- 
erage attrition length of kinetic growth SAWs, at which 
they cannot continue without violating the self-avoidance 
condition. The paper closes with the conclusions in 
Sec. V. In three appendices we give details of the proba- 
bilistic calculations. 



II. DEFINITIONS AND MODELS 



properties with system size, we simulated networks of 
several sizes, including up to 3500 x 3500 sites. Periodic 
boundary conditions were assumed. 

For small-world networks, the probability distribution 
of connectivities, P sw (k), is short-tailed with an expo- 
nential decrease for large k. Analytical expressions for 
this probability distribution have been given elsewhere 
[201 1431] . and will not be repeated here. We will only give 
some details necessary for the discussion below. To this 
end, we distinguish now between the ends of links that 
remain on their original sites (as in the regular lattice), 
and those changed in the rewiring process. The number 
of links in a network with mean connectivity z amounts 
to zN/2. Since each link connecting two sites in the 
starting regular lattice is rewired with probability p, the 
average number of changes per site is zp/2. This means 
that each connection of a site is removed in the rewiring 
process with probability p/2. Then, the number of links 
s generated in the rewiring process and arriving at a node 
follows, in the limit of large N, a Poisson distribution for 
s > 0: 



P r (s) = 



= 1 (wy p - v /2 



(i) 



We note that other ways of generating small- world net- 
works from regular lattices have appeared in the litera- 
ture. In particular, Newman et al. [19l l44j instead of 
rewiring each bond with probability p, added shortcuts 
between pairs of sites chosen uniformly at random, with- 
out removing any bonds from the starting lattice. This 
procedure turns out to be more convenient for analytical 
calculations, but does not keep constant the mean degree 
(k) , which in this case increases with p. 



A. Small-world model 

The networks studied here were generated from a two- 
dimensional square lattice (connectivity z = 4). Small- 
world networks were built up according to the model of 
Watts and Strogatz [y, [lQ|. We consider in turn each 
of the bonds in the starting lattice and replace it with 
a given probability p by a new bond. This means that 
one end of the bond is changed to a new node chosen at 
random in the whole network. We impose the conditions: 
(i) no two nodes can have more than one bond connecting 
them, and (ii) no node can be connected by a link to itself. 
This method keeps constant the total number of links in 
the rewired networks, and the average connectivity (or 
degree) is (fc) = z. The total number of rewired links is 
zpN/2 on average. In the rewiring process we avoided 
sites with zero and one links, and thus each site has at 
least two neighbors. Note that nodes with k = 1 behave 
as culs-de-sac for self-avoiding walks, i.e., a walk arriving 
at a node with connectivity k = 1 cannot continue, even 
though it has not yet intersected itself. 

To study the asymptotic behavior of kinetic growth 
SAWs as N — > oo, and analyze the convergence of their 



B. Self-avoiding walks 

A self-avoiding walk is defined as a walk along the 
bonds of a given network which can never intersect itself. 
The walk is restricted to moving to a nearest-neighbor 
site during each step, and the self-avoiding condition con- 
strains the walk to occupy only sites which have not been 
previously visited in the same walk. The simplest pro- 
cedure to obtain SAWs consists in generating ordinary 
random walks and stop when they arrive at a node al- 
ready visited. In regular lattices, there appears a problem 
with this sampling algorithm, caused by the rapid attri- 
tion for long walks, because the probability of an n-step 
walk being self-avoiding decays exponentially for large 
n [45j. Due to this important limitation, more complex 
methods based on Monte Carlo sampling have been used 
to generate SAWs with the correct weight, and to obtain 
ensemble averages of various quantities [45| . In general, 
for networks including nodes with different degrees (con- 
trary to usual regular lattices), sampling by employing 
simple random walks introduces a bias in the weight of 
different SAWs. 




FIG. 1: Schematic diagram showing a non-reversal random 
walk of length n = 5 on a network. Open and black cir- 
cles represent unvisited and visited nodes, respectively. The 
starting node is indicated as 0. After arriving at site 5, the 
non-reversal condition allows to choose for the next (sixth) 
step one of three possible nodes (labeled as 1, 6, and 7). For 
a non-reversal SAW, one chooses among these three nodes. If 
1 is selected, then the walk stops. For a kinetic growth walk, 
one chooses 6 or 7, each with 50% probability. 



walks [45j . In these walks one randomly chooses the next 
step among the neighboring nodes, excluding the pre- 
vious one. If one chooses a node already visited, then 
the walk stops (see Fig. 1). These walks will allow us 
to study the self-intersection length (see Sec. III). The 
second type of walks considered here are kinetic growth 
walks [49], in which one randomly selects the next step 
among the unvisited adjacent sites and stops growing 
when none are available. These walks have been stud- 
ied earlier to describe the irreversible growth of linear 
polymers [49j . and will allow us to consider the attrition 
length for a walk on a given network (see Sec. IV). We 
note that kinetic growth walks are less sensitive to at- 
trition than non-reversal SAWs, since in the former the 
walker always escapes whenever a way exists. 

For each set of parameters (p, N) , we considered differ- 
ent network realizations, and for a given network we se- 
lected at random the starting nodes for the SAWs. For a 
given parameter set, the total number of generated SAWs 
amounted to 10 6 for each type of walks. 



III. SELF-INTERSECTION LENGTH 



For regular lattices, the number u n of SAWs starting 
from a site has an asymptotic dependence for large n 
[IE El: 



/' 



(2) 



where 7 is a critical exponent which depends on the lat- 
tice dimension, and /i is the connective constant or effec- 
tive connectivity of the considered lattice [4a,|43]. In gen- 
eral, for a lattice with connectivity z, one has // < z — 1. 
This parameter fi can be obtained from Eq. ([2]) as the 
limit fi = linin^oo u n /u n -\. The connective constant de- 
pends upon the particular topology of each lattice, and is 
known with hig h accuracy for two- and three-dimensional 
lattices [45|, [48( . Self-avoiding walks on small- world net- 
works have been studied in recent years [381 ]. In partic- 
ular, for networks generated from the square lattice, the 
connective constant fi was found to rise from fi = 2.64 for 
rewiring probability p — (regular lattice) to /i = 3.70 
for p = I. 

To analyze dynamical processes on networks, it can be 
more convenient to use kinetically grown SAWs. These 
walks are well-suited to study search or navigation pro- 
cesses on networks, where they are assumed to grow step 
by step in a temporal sequence. Their asymptotic be- 
havior can give us also information about the long-range 
properties of complex networks, in particular on the pres- 
ence of loops. In conection with this, an interesting prop- 
erty of kinetically grown SAWs is the number of steps 
necessary to intersect themselves. Moreover, we are in- 
terested in the possibility of avoiding visited nodes and 
maximize the length of the walks. For these reasons, 
we will consider in the following two kinds of growing 
walks. The first kind will be non-reversal self-avoiding 



To study the probability of a walk intersecting itself, 
we consider non-reversal walks that stop when they reach 
a node already visited in the same walk. The number of 
steps of a given walk before intersecting itself will be 
called self-intersection length and will be denoted I. 

Let us consider Mq non-reversal walks starting from 
nodes taken at random, and call M(n) the number of 
walks surviving after n steps (i.e., those which did not 
arrive yet at any node visited earlier). In regular lat- 
tices, it is known that the ratio M(n)/Mo decreases with 
the number of steps n as n 7_1 e _A °™, where 7 is the same 
exponent as in Eq. @, and Ao = log[(z — l)//x]. Hence, 
the surviving probability basically decays for large n as 
exp(— Xqu), with Ao = 0.129 for the two-dimensional 
square lattice [451 ]. 

We now turn to the results of our simulations for small- 
world networks. In Fig. 2 we show the fraction of remain- 
ing walks, M(n)/Mo, as a function of n for networks with 
several sizes and rewiring probability p = 0.3. The net- 
work size increases from left to right, and one clearly 
observes a finite-size effect on the surviving probability 
of the walks. For increasing network size, the results con- 
verge to an exponential decrease of the ratio M(n)/Mo. 
In fact, for the maximum size presented in Fig. 2 
(2000x2000 nodes), we find M(n)/M = exp(-An), with 
A = 0.014. For other p values, one finds similar decays 
of the number of surviving walks, with the coefficient A 
decreasing as p rises. 

Let us call t n the conditional probability of a walk 
stopping at step n+1, assuming that it actually arrived 
at step n. The exponential decay observed for M(n)/Mo 
is compatible with a constant t n , independent of n (i.e., 
t n = £), as is expected at least for the asymptotic regime 




200 400 600 800 

Number of steps n 

FIG. 2: Fraction of non-reversal SAWs that survive after n 
steps, without intersecting themselves. Results are plotted for 
small-world networks with rewiring probability p = 0.3 and 
different sizes, as derived from numerical simulations. From 
left to right: L = 100, 200, 300, 500, and 2000. 



at large n. In fact, assuming a constant t n , one has 
M(n) - M(n + 1) = tM(n) , 



(3) 



which can be solved by iteration with the initial condition 
M(0) = M , giving 



M(n) = M (l - t) r 



(4) 



Comparing this result with the exponential decay of 
M(n)/M found above, we can identify: A = — log(l — t). 
The mean self- intersection length {l) sw in small- world 
networks is plotted in Fig. 3 as a function of the system 
size N for several values of p. For a given N, (l) sw in- 
creases as p is raised, as could be expected. For a given 
rewiring probability p, (1} SW rises with the system size, 
and eventually saturates to a finite value for large N. 
This saturation appears at larger system sizes as p is in- 
creased. Finally, for p = 1 we find a behavior similar 
to that expected for random networks, for which (l) rn 
increases as a power of N [4fJ : 



(D 2 « — 
[)rn 2w 



(5) 



where wis a network-dependent parameter, w = ((k 2 ) — 
2(k))/(k) 2 . The dashed line in Fig. 3 corresponds to 
Eq. (JSJ) for a random network with minimum degree 
fc = 2 and (k) = 4, which gives (l) rn = 1.597V 1 / 2 [4(J- 
Thus, the results of our simulations for small-world net- 
works with p — 1 agree with those expected for random 
networks, for which {l) sw diverges with the system size 
asiV 1 / 2 . 
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FIG. 3: Average self-intersection length {l) sw as a function of 
system size N for small-world networks with several rewiring 
probabilities p. Symbols display results of numerical simula- 
tions: open squares, p = (regular lattice); circles, p = 0.2; 
triangles, p = 0.4; diamonds, p — 0.6; filled squares, p — 1. 
Error bars are less than the symbol size. Dotted lines are 
guides to the eye. The dashed line shows an analytical pre- 
diction for random networks, given by Eq. ([5]). 



The mean self-intersection length in the thermody- 
namic limit is presented in Fig. 4, as derived from the re- 
sults of our simulations (open symbols). In the remainder 
of this section, we present two analytical approximations 
for (l)sw, valid for the regions close to p — and p = 1. 
These approximations correspond to the lines shown in 
Fig. 4. 

To derive an approximate expression for the mean self- 
intersection length for finite (small) p, we consider non- 
reversal SAWs on the regular lattice, and calculate the 
corrections introduced by the rewired links in small-world 
networks. An important parameter in this calculation is 
the conditional probability q of following a rewired link 
in step n + 1, assuming that in fact the walk reached step 
n. To first order in p, this probability is given by (see 
Appendix A): 



P 



z-1 



(6) 



where the term z — 1 refers to the number of available 
connections for a given step n of a non-reversal walk on 
the regular lattice. For the derivation presented in Ap- 
pendix A, we assume: (1) The walks include zero or one 
rewired links, which is compatible with our order-p ap- 
proximation, and (2) a rewired link takes the walker far 
away from the starting node, which is valid in average 
for large system size. 

With these assumptions, we obtain for small p 
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Rewiring probability p 

FIG. 4: Mean self-intersection length {l) S w as a function 
the rewiring probability p in the thermodynamic limit N — > 
oo. Symbols represent results of numerical simulations. Error 
bars are less than the symbol size. Lines indicate results of 
analytical calculations for p <^ 1 (dotted line) and for p near 
1 (dashed line). 



Eq. 



(l)s W »(l)2D + U(l)2D + 2(l) 2 2D 



(l 2 hn) 



(7) 



where the subscript 2D indicates that the average values 
are taken for walks on the square lattice. Thus, the mean 
self-intersection length of non-reversal walks on small- 
world networks with p -C 1 is expressed in terms of the 
rewiring probability p (related with q by Eq. ([6])) and 
known results of the regular lattice. The mean length 
(l) sw given by Eq. ([7]) is plotted in Fig. 4 as a dotted 
line. 

We now consider the average self-intersection length 
for rewiring probability p near 1. As shown above, at 
p = 1 we find for small-world networks the behavior typ- 
ical of random networks, with (l) rn increasing as A 1 / 2 , 
and eventually diverging in the thermodynamic limit. 
For small- world networks with p < 1, {l) sw has a finite 
limit for N — > oo, as shown in Fig. 3, and its value is de- 
termined by correlations present in the networks. Even 
though these correlations are small for p close to 1, they 
still control the long-range behavior of SAWs. In fact, 
in the limit p — > 1 and for large enough N, the factor 
limiting the length of the walks is the residual presence 
of loops of size four in the networks, remaining from the 
starting regular lattice. Thus, we assume that the walks 
finish after circulating along the links of four-node loops 
(see Fig. 5(a)), where no one of the connections forming 
them was rewired. Taking into account that in the square 
lattice, the number of loops of size 4 equals the number 
of nodes A, then the average number N$ of remaining 
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FIG. 5: Schematic representation of SAWs along four-node 
loops in small-world networks, (a) A non-reversal walk ar- 
riving at site 1 and following through sites 2, 3, and 4. The 
walk stops if it tries to return to site 1 (dashed arrow), (b) 
A kinetic-growth walk arriving at a four-node loop including 
a site with k = 2 (labeled as 4). The walk stops if it follows 
the sequence of nodes 1, 2, 3, 4, since at node 4 it cannot 
continue. 



four-node loops is 



As 



N(l~p) 4 



(») 



Note that in the rewiring process new loops of size 4 (and 
also of size 3) can appear, but for large N, their number is 
much less than As. In fact, the average number of four- 
node loops in random networks is basically independent 
of A: A 4 w z(z-l) 3 /2, and thus for any p < 1, N s > A 4 
for large enough A. 

With these assumptions, we can calculate the condi- 
tional probability t of a non-reversal SAW stopping at 
step n+1, assuming that it arrived at step n. This is done 
in Appendix B. The probability t turns out to be inde- 
pendent of n for large system size, and then for p near 1, 
the probability distribution Q{1) for the self-intersection 
length I can be approximated as 



Q(l) = t(l-t) 1 



(9) 



which gives the probability for a walk proceeding until 
step I, and being stopped at step I + 1. 

From the distribution Q(l), one finds an average self- 
intersection length (l) sw = (1 — t)/t. For p close to 1, one 
has t <C 1, so that, using the expressuion for t derived in 
Appendix B, we find 



{'), 



81 



(1-P)- 4 . 



(10) 



This result takes only into account the self-intersection 
at loops of size 4, which gives a dependence of t in the 
form (1 — p) 4 . Additional higher-order terms of the form 
(1 — p) m can be obtained by considering loops of size m 
in the analyzed networks (m — 6, 8, ...), but their con- 
tribution will be negligible close to p — 1 . The analytical 
dependence of (l) sw given by Eq. (|10[) is shown in Fig. 4 
by a dashed line. It shows good agreement with results 
of the numerical simulations for p > 0.5, and predicts a 
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FIG. 6: Mean attrition length {L) sw as a function of system 
size N for small- world networks with several rewiring proba- 
bilities p. Symbols represent results of numerical simulations: 
circles, p = 0.2; triangles, p = 0.4; diamonds, p — 0.6; filled 
squares, p — 1. Error bars are smaller than the symbol size. 
Dotted lines are guides to the eye. The dashed line is an ap- 
proximate analytical prediction for random networks, given 
in Ref. [4Q] (see text for details). 



divergence of {1} SW for p — > 1, in line with the behavior 
known for random networks. On the other side, for any 
p < 1 our calculations predict a finite value for (1} SW . 

It is worth commenting that using the same approx- 
imations employed above, one finds that the standard 
deviation of the self-intersection length of non-reversal 
SAWs is <Ji rs 1/t, or ai » (l) sw . This is consistent with 
the exponential behavior of the surviving probability of 
these walks, as discussed above. 



IV. ATTRITION LENGTH 

In this section we consider random SAWs, that grow 
on the network until they arrive at a node (called blocking 
node in the sequel), where they cannot continue because 
all adjacent nodes have been already visited and are not 
available for the walk. These are kinetic growth walks, 
as introduced by Majid et al. [4!|, and defined above in 
Sec. II. The number of steps of a given walk until being 
blocked will be called attrition length, and will be denoted 
L. 

In Fig. 6 we show the mean attrition length (L) sw as 
a function of the system size N for several values of the 
rewiring probability p. Symbols represent results of nu- 
merical simulations, and dotted lines are guides to the 
eye. For a given N, (L) sw increases as p rises, simi- 
larly to the mean self- intersection length (see Fig. 3). 
For a given rewiring probability p < 1, (L) sw rises with 
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FIG. 7: Mean attrition length {L) sw as a function of the 
rewiring probability p in the thermodynamic limit N — ► oo. 
Symbols represent results of numerical simulations. Error 
bars are less than the symbol size. Lines show results of an- 
alytical calculations in the small-p limit (dotted line) and for 
p close to 1 (dashed line). 



the system size, and saturates to a finite value for large 
N . For p = 1, however, we obtain an increase of (L) sw 
as a power of N in the whole range of network sizes 
studied here (up to N s=s 10 7 sites). In fact, we find 
(L) sw ~ N a , with an exponent a = 0.51 ± 0.01. The 
dashed line in Fig. 6 represents an analytical approxima- 
tion for the mean attrition length in random networks (as 
given by Eqs. (10) and (24) in Ref. 40]). According to 
this calculation, the mean attrition length in uncorrelated 
networks with short-tailed degree distribution scales as 
(L) rn ~ iV 1_1 / fc °, where fco is the minimum degree in the 
network. For the networks considered here, ko = 2, and 
one should expect (L) sw ~ A 1 ' 2 , in line with the results 
of our simulations shown in Fig. 6 (filled squares). 

We note that the small-world networks with p = 1 
studied here are not totally uncorrelated, because the re- 
striction fco = 2 imposes that some correlations coming 
from the starting regular lattice are still present in the 
rewired networks. These correlations seem to be unim- 
portant for defining the main trends of the long-range 
behavior of kinetic-growth walks. 

Shown in Fig. 7 is the mean attrition length in the 
thermodynamic limit, as derived from our simulations 
of kinetic-growth SAWs on rewired networks (open sym- 
bols). For the system sizes considered here, it is hard to 
obtain a definite value for the limt N — ► oo for rewiring 
probabilties p larger than 0.6. Lines in Fig. 7 indicate 
results of analytical approximations presented below. 

To calculate the mean attrition length of kinetic- 
growth walks on small-world networks with p close to 
cero, we employ a method similar to that used for the self- 



intersection length in Sec. III. In fact, one can follow ex- 
actly the same steps as in Appendix A, but putting only 
the connective constant ^ instead of z — 1 in Eqs. (jAlj) 
and (|A2[) . In the non- reversal SAWs employed to cal- 
culate the self-intersection length, one can choose among 
z — 1 sites at a given step n on the regular lattice. Instead 
of this, for kinetic-growth walks used to obtain the attri- 
tion length, one has in average /i sites available for step n 
(assuming n>l). For the square lattice employed here 
to obtain the small-world networks, we have z — 1 = 3 
vs /i = 2.64. Then, we find for the mean attrition length 
for small p: 



(L) sw = (L) 2D + | [(L) 
where 



>o + 2(L) 2 D -(L 2 ) 2D ] 



(11) 



(12) 



The dependence of (L) sw on p given by Eqs. (fTTj) and 
(fT2f is shown in Fig. 7 as a dotted line. As expected, 
this approximation to order p gives values of (L) sw lower 
than the numerical simulations when p is increased from 
p = 0. This is due to the fact that it neglects higher-order 
terms, with the contribution of walks going on two, three, 
..., rewired links. This contribution can only increase the 
mean attrition length given in Eq. pip , which takes into 
account walks going over zero or one rewired connection. 

To obtain the asymptotic behavior of the mean attri- 
tion length (L) sw close to p — 1, we follow a procedure 
analogous to that employed to calculate (l) S w Here, as 
above, {L) sw is limited by the presence of loops in the 
network, and in particular by loops of size four remain- 
ing from the starting square lattice. This means that 
the highest probability for a kinetic-growth walk being 
blocked occurs at four-node loops, but in this case there 
appears another condition for a path to be stopped, since 
one needs a blocking node with degree k = 2 (if all four 
nodes in a loop have k > 2 the walker will escape, see 
Fig. 5(b)). 

With these assumptions, we can calculate the (con- 
ditional) probability t' of a walk being blocked at step 
n + 1, assuming that it actually reached step n. For 
large n, t' can be considered independent of n and then 
the distribution of attrition lengths can be expressed as 
Q'(L) = t'(l — t') L . Finally, the mean attrition length 
can be approximated for f < 1 as (L) sw «s 1/t', and 
using the expression for t' derived in Appendix C, we 
find 



(L) sw ^27e 2 Pp- 2 (l- P y 



(13) 



This dependence of {L) sw on the rewiring probability p is 
displayed in Fig. 7 as a dashed line. The trend predicted 
by Eq. (fl~3"|) is similar to that of the results of numerical 
simulations in the region 0.4 < p < 0.6, but somewhat 
higher (less than a factor of 2). Close to p = 1 we find a 
divergence of (L) sw as (1 — p)~ 4 , like in the case of the 
self-intersection length (l) sw (see Eq. (fTU)) ). 



For small-world networks rewired from a cubic lat- 
tice, we expect for p — » 1 divergences of the mean self- 
intersection and attrition length similar to those found 
here for networks derived from a two-dimensional square 
lattice. In fact, for the cubic lattice the minimum-size 
loops include four nodes, and those remaining after the 
rewiring process will control (1} SW and {L) sw close to 
p = 1, which will diverge as (1 — p)~ 4 . Something anal- 
ogous will happen for networks rewired from hypercubic 
lattices of dimensions higher than three. 



V. CONCLUSIONS 

Self-avoiding walks provide us with an adequate tool to 
analyze long-range characteristics of complex networks. 
In particular, they allow us to study the quality of a net- 
work to be explored without returning to sites already 
visited. Here, we have studied the self-intersection length 
I and attrition length L of kinetically-grown SAWs on 
small- world networks, rewired from a two-dimensional 
square lattice. With this purpose, we have considered 
non-reversal SAWs to obtain (l) sw and kinetic-growth 
walks to obtain (L) sw . 

We have calculated self-intersection and attrition 
lengths by means of approximate probabilistic methods, 
which give results in good agreement with those derived 
from numerical simulations. For rewiring probability 
p = 1, both the average self- intersection and attrition 
length diverge with the system size as N 1 / 2 . This de- 
pendence of (L) sw , however, changes with the minimum 
degree fco present in the networks, and in general (L) sw 
scales as N a with an exponent a = 1 — 1/fco (as for ran- 
dom networks). 

In the thermodynamic limit, both {l) sw and {L) sw in- 
crease with p, but remain finite for p < 1. The length 
of kinetically-grown SAWs is limited by the presence of 
loops in the networks, mainly those loops remaining from 
the starting regular lattice. Close to p — 1, both mean 
lengths diverge as (1 — p)~ A . The same type of divergence 
is expected to appear for small-world networks rewired 
from cubic and hypercubic lattices in higher dimensions. 
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APPENDIX A: MEAN SELF-INTERSECTION 
LENGTH FOR p -> 

Here we calculate the mean self-intersection length for 
small p from the known results for p = (regular lattice). 
To this end, let us consider a non-reversal SAW on the 



8 



regular lattice, that has reached at least length n. Let us 
call S n the node visited in step n. For step n + 1 on the 
small-world network, we will distinguish three different 
possibilities, depending on the nature of the link followed 
in this step: 

(1) The link employed in step n + 1 was rewired, but 
keeping one end on node S n (probability p/2). 

(2) The path follows a link that was rewired from a 
distant node to node S n . Since each node is on average 
connected in the rewiring process to pz/2 = 2p new links, 
the conditional probability of following one of these links 
in step n + 1 is 



P 



2p 



z-l + 2p 



(Al) 



which gives to order p: p' ss 2p/(z— 1). In this expression, 
z — 1 is the number of available links for a step of a non- 
reversal walk on the regular lattice. 

(3) The link followed in step n + 1 was not rewired, and 
remains as in the starting regular lattice. This occurs 
with probability 1 — q, where q — p/2 -\-p' , or to order p: 



P 



z-\ 



(A2) 



Then, q is the probability of following a rewired link in 
step n+1, assuming that in fact the walk reached step 
n. 

Now, let us consider a non-reversal SAW of length I on 
the regular lattice, and assume that it intersected itself 
at step I + 1. The same walk will still be possible on a 
rewired network (when no one of the links it uses was 
rewired) with probability (1 — qy. For a given n < I, 
the average length of walks going on a rewired connec- 
tion at step n + 1 is n + {I) 2D, where {V)2D is the mean 
self-intersection length in the regular lattice. This in- 
cludes the assumption that a rewired link takes in gen- 
eral the walker to a node on the lattice far away from 
the starting one (which is valid in average for large N). 
We also assume that once a walk has employed a rewired 
link, it does not find any other, which is consistent with 
our order-p approximation for (l) sw (no more than one 
rewired link in a walk). Therefore, for walks of length / 
on the regular lattice, the average length (l')i for walks 
on the rewired networks is: 



(l') l = l(l-q) l + J2^+(l)2D)qn 



(A3) 



where the second term on the r.h.s. takes into account 
the contribution of walks going on a rewired link, and the 
first, that of walks visiting no rewired links. Here q n is 
the probability of a walk going on a rewired link at step 
n, given by q n = q(l - q) 71 ' 1 . 

Then, to first order in q we have: 



{l')l~l + ql 



l-l 



(I) 



2D 



(A4) 



Finally, we calculate the average of (/'}/ over all possible 
walk lengths I in the regular lattice, and find: 



(I) aw ~ (l)iD + Q 



~ 2 dhD 



(1)Id 



\(?hD 



(A5) 



APPENDIX B: MEAN SELF-INTERSECTION 
LENGTH FOR p -► 1 

Here we present an approximate calculation for {l) sw in 
small- world networks with p close to 1. In this limit, and 
for large enough N, the length of the walks is basically 
limited by the residual presence of loops of size four in 
the networks. The number of these loops remaining from 
the starting square lattice is given by Eq. ([5]): N$ = 
N(l-p) 4 . 

We then assume that the walks finish after circulating 
along the links of four-node loops (see Fig. 5(a)). Given 
a loop of size 4 in a small-world network, for each of its 
nodes there are in average two (directed) links leading 
to the loop from the outside. Then, the total number of 
links leading to nodes in these loops is 



A ; =8A(l-p) 4 



(Bl) 



Since the number of directed links in the network is AN, 
then the fraction of those leading to nodes in loops of 
length 4 is / = 2(1 — p) 4 . Note that we consider nondi- 
rected networks, but for our present purpose we have to 
distinghish between the two ways in which a single link 
can be visited in a walk. 

Let us now assume that a walk has arrived at a node 
(say the node labeled "1" in Fig. 5(a)) in a four-node loop 
in step n — 4 (n > 4). In this case we look for the con- 
ditional probability of circulating around the loop and 
returning to node "1" after four steps, i.e., in step n. 
This probability r is given by r = 2(l/3) 4 = 2/81, where 
we have used that the mean connectivity is (k) = 4, 
and that there are two ways to walk around the con- 
sidered four-node loop (clockwise and counter-clockwise, 
see Fig. 5(a)). Finally, the conditional probability t of 
arriving in step n at a node visited earlier (in step n — 4) 



t = rf = — (1-p) 4 
.1 gl V V) 



(B2) 



Note that t defined in Eq. (|F32|) coincides for p close to 1 
with the stopping probability introduced in Eq. Q, and 
is independent of the step n (for n > 4) . 



APPENDIX C: MEAN ATTRITION LENGTH 

FOR p -> 1 

In this appendix we calculate the probability of ar- 
riving at a node with k — 2 (labeled as "4" in Fig 
5(b)) in a four-node loop, after circulating along the 



links of the loop. We assume that node "4" has de- 
gree two, so that two of its original links in the regu- 
lar lattice have been rewired away (probability p/2 for 
each one), and no other connection has been established 
with this node during the rewiring process (probability 
P r (0) = cxp(— 2p), see Eq. (p}). Hence, the average num- 
ber of nodes with k = 2 in the remaining four-node cycles 
is N v = AN s {p/2) 2 P r (0), with N s given in Eq. ©. This 
yields 



N v = JV(1 - p) V e 



-2p 



(CI) 



In order for a walk to stop at node "4" (in Fig. 5(b)), 
it has to enter into the loop through nodes "1" or "3", 
and then circulate along the loop counter-clockwise or 
clockwise, respectively. Since the average number of links 
leading to the loop from the outside through nodes "1" 
or "3" is four, the fraction of links in the network leading 
to this kind of cycles is 

4N, 



/'=^ = (l-P)Ve 



-2p 



AN 



(C2) 



This includes the assumption that the probability of four- 
node cycles with more than one node with k — 2 is negli- 
gible (in fact it is much lower than that of only one node 
with degree 2). 

If a walk has arrived at nodes "1" or "3" from the out- 
side of the loop, then the conditional probability of cir- 
culating around the loop and arriving at node "4" (the 



blocking node) is 



(1/3) , i.e. 1/3 for each of three 



successive steps along the loop. Thus, similarly to Sec. 
Ill and Appendix B, the conditional probability if of ar- 
riving in step n at a node with k = 2, after circulating 
around a four-node loop (in the form shown in Fig. 5(b)), 
is t' = r'f, or 



t' = ^e-^{l-pf 



(C3) 
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